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We study the effect of a background flux string on the vacuum energy of massive Dirac fermions 
in 3+1 dimensions confined to a finite spatial region through MIT boundary conditions. We treat 
two admissible self-adjoint extensions of the Hamiltonian. The external sector is also studied and 
unambiguous results for the Casimir energy of massive fermions in the whole space are obtained. 
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(S| , I- INTRODUCTION 

The influence of background fields or of boundaries on the vacuum structure of a quantum field theory is of 
continuing importance in various branches of modern physics. The various aspects are often described generically as 
Casimir effects and many techniques have been developed in order to analyse the different situations. As examples 
■ we mention the Green's function approach Jl|-|j| and the zeta function regularization [§-0 . 

In the years since the Casimir effect was first discussed Q these techniques have been refined considerably. Relatively 
recently, a contour integral method has been developed which allows the representation of the Casimir or ground state 
energy in terms of eigenfunctions or of the Jost function of the associated field equation |p[-pj|. This allows for a 
relatively direct analysis for configurations where this information can be easily obtained. As a result, spherically or 
cylindrically symmetric situations have been analysed in various contexts, see e.g. pd|-p^|. Most of the research done 
so far concentrates on the influence of background fields or of boundaries separately. Relatively little is known when 
these two effects are combined. Exceptions are e.g. [ p0| , ^6[ , where a magnetic fluxon and MIT boundary conditions 
t-H ' for a massless and massive Dirac field in 2 + 1 dimensions were considered. 

It is our aim to continue the analysis of this combined effect in a 3 + 1 dimensional cylindrically symmetric situation. 
We consider a massive Dirac operator in the presence of a magnetic flux string located along the z-axis and with MIT 
boundary conditions at a cylinder of radius R. Because of the presence of the flux string, a self-adjoint extension 
of the radial Dirac operator is needed, as it is well known from the 2 + 1 dimensional situation |^l],^2) . In analogy 
to (ilf | we will consider two possible self-adjoint extensions which are both compatible with the presence of a Dirac 
1 delta magnetic field at the origin. One extension is constructed by imposing spectral boundary conditions ]23]-|25|] at 
a finite radius which is then shrunk to zero. The second one follows from the zero radius limit of a cylindrical flux 
shell ]26|-|28|] . The Casimir energies differ, reflecting the fact that different self-adjoint extensions describe different 
physics in the core. 

Compared to the calculations in 2 + 1 dimensions, peculiar differences occur. Most importantly, when considering 
only the interior of the cylinder, the poles in the Casimir energy depend on the flux for both extensions. This renders 
an interpretation of the results impossible, in the sense that finite values cannot be extracted in a physically reasonable 
way. This led us to consider also the exterior space, in an attempt to get a simpler pole structure by adding up both 
contributions. In doing so, unambiguous results for Casimir energies in different situations can be obtained. Let us 
stress that in the present situation general heat-kernel arguments [£9|,[50| can not be applied simply because no general 
answers are known for cases with singular fluxes and boundaries. For this reason we really need to perform explicit 
calculations. 

The paper is organized as follows. First we briefly present the solutions of the Dirac equation in the presence of a 
flux string. In the following sections self-adjoint extensions are obtained and some details of how to impose spectral 
boundary conditions are provided. Having chosen the self-adjoint extensions, the MIT boundary conditions are 
imposed and implicit eigenvalue equations are obtained. These are the basis for the analysis of the Casimir energy. 
In order to obtain a well-defined Casimir energy (such that a numerical analysis makes sense) we have evaluated 
differences between the case when the flux is arbitrary and when the flux is integral. Also, we include the exterior 
space and the pertinent changes in the calculation are briefly explained. A numerical analysis of Casimir energies as 
a function of the flux, for the Dirac field, is presented afterwards. 
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II. DESCRIPTION OF THE PROBLEM 



We study the Dirac equation for a massive particle in 3+1 dimensional Minkowski space-time in the presence of a 
flux string located at the origin, 



H = V A A = -8{r)e z 
r 



(1) 



We are going to impose boundary conditions on a cylinder and for that reason use cylindrical coordinates. The gamma 
matrices in the chiral representation are 

J° = P3®V3, 7 1 = «P3 ® cr 2 , ; 7 2 = -ip 3 (g) CTi, ; = ip 2 ® I2, 
and in a comoving coordinate frame one has the usual relation 
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cos 6 j 1 + sin6>7 2 , ~f = - (-sin #7* + cos#7 2 ) 



Furthermore, for the flux 
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— , k = I + a, < a < 1. 
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Here we have introduced the integer part I of the reduced flux k such that a is restricted to the given interval. 
Our main concern is the Casimir energy and for this reason we consider the Dirac Hamiltonian for the system, 



H B ^E = -*7V [~ d r + Yr(d e - in) + 'f'fd, + imf] * 



(2) 



The representation of the Hamiltonian is chosen in such in a way as to simplify the implementation of spectral 
boundary conditions at a later stage. Given the cylindrical symmetry of the configuration the eigenfunctions have 
the form 



— ik z z 



For ijje{r, 0) there are four types of solutions. They are 
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where n = — 00 , 00, and we have used the notation 

k = \/ e 1 — to 2 , 



S = sign(e)\A 2 + e 2 



(3) 



(4) 
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-— (e — -E 1 ) , w = n — I — a. 

k z 



Here e 2 = A: 2 + m 2 are the eigenvalues of the square of the two-dimensional Dirac operator obtained after setting 
k z = 0. 



In order to fix the eigenvalues E we have to choose a self-adjoint extension of the radial Hamiltonian and we have 
to impose boundary conditions. A family of self-adjoint extensions arises and we will consider two particular cases. 
The first one is obtained by imposing spectral boundary conditions at some interior cylinder the radius of which is 
sent to zero. Spectral boundary conditions are imposed as described in p3]-p5[|, see also [ pi] ]. In the present case a 
suitable choice for the boundary operator A is [pl|,p2[ 
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(6) 



The term (l/2r)l4 = (l/2r)Kl4, with K the extrinsic curvature of the boundary, has been added to guarantee that 
the operator A leads to a self-adjoint boundary value problem, see |n| for details. 

Spectral boundary conditions amount to setting to zero the projection of ^ e onto all cigcnfunctions of A with 
negative eigenvalues at the boundary. Again, due to the symmetry, the eigenfunctions have the form 



A = 



—ik z 



(7) 



where again four different types of solutions exist, 
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Here we have introduced 
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and «i, «3 can be seen to have eigenvalues 0.1,3 = 0-/ r o) a + an d 0^2, «4 have instead 111,4 = ~(l/ r o)a-- The 
eigenvectors and eigenvalues of A in 2 + 1-dimensions are obtained by taking fc 2 — > 0. 

Having these eigenfunctions a at our disposal, projections with can easily be performed. Due care has to be 
taken in order to implement the boundary conditions correctly for all values of n. In the limit ro — > 0, the vanishing 
of the projection rules out some of the ^4 in Eqs. (0) and (@) because they diverge at the origin. In detail, one 
finds that for n > I + 1 the function ipi™' 1 ^ and ipi™' 3 ^ form a suitable basis, whereas for n < I the relevant ones are 
The remaining one, n = /, allows the identification as a member of the one parameter family of 



-0e™' 2 ^ and tp, 
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self-adjoint extensions 



Writing the spinors in the form 
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self-adjoint extensions are classified according to the conditions 
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In our case, for < a < 1/2, 
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and, for 1/2 < a < f, 



/rw = 
srM = 



[-A+C- +rAfC+] J- a+1 (kr ), 
[A+r + A-] J_ Q (fcr ), 
[A+rC_- J 4 ; -C+] J_„+i(Ato), 



[B++TB-] J a (kr ), 
[B+C--TB-C+] J a -i(kr ), 
[B+T + B-] J a {kr ), 
[-B+TC-+B-C+] J a -i(kr ) 
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Considering the above condition (^) this imposes 



e ± = l + f ° <a< f . (io) 

We refer to this choice as I. Our second choice, called II, corresponds to 

This extension has been considered in [p7f and it arises when a finite radial flux, which is ultimately shrunk to zero, 
is considered. 

We are now in a position to impose MIT boundary conditions at the exterior boundary r = R. These boundary 
conditions ensure that the fermion current across the boundary vanishes. The relevant boundary operator is 

; ( J o I )■ <'-> 



with 



and the boundary conditions reads 



B^ E {R,9,z) = 0. 



Given the previous discussion for type I, see Eqs. (||) and (|i"o|), this boundary condition has to be applied on the 
following set of spinors, 



n> l + l : i>^ l \ Vi"' 3) , 
n<l-l; Vi"' 2) , 4 nA) , 
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and similarly for type II. As a result, the eigenvalues take the form 

1 



E=-^(k z R) 2 +z 2 + x 2 (15) 
H 

where z = mR and x are the solutions of the equations obtained from 

- - ^Mx)J^-i(x) - 0, (16) 

with /j, = v ± a, v = n + 1/2 = 3/2, 5/2, . . . , oo and a = a — 1/2. The critical subspace v — 1/2 belongs to the case 
[i = v — a for— 1/2 < a < and to the case fi = v + a when < a < 1/2. All eigenvalues have degeneracy two 
coming from particle and antiparticle states. 

III. THE CASIMIR ENERGY FOR THE INTERIOR 

Now everything is prepared for the calculation of the Casimir energy in the zeta function regularization scheme. 
Our presentation shows only the Type I case. We comment on the the minor changes needed for Type II where 
appropriate. Given the translational invariance in the z-direction, we define the Casimir energy density by 

4^) = _^M 2 Xi„t( S )| s= -i/ 2 , (17) 
where the zeta function has the structure 
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Expressed as a sum over the zeros of Eq. (|16|) , we can write 

Cint( S )=ECf t) ( S ) ( 19 ) 



with the partial zeta functions 



c m )(s) = ^ 2s -ir^_V2) g (z2 + x , im) -(-i/2) (20) 

v ^ * rn—1 

Slightly more explicit is the form 

oo 

Cint( S ) = -C$!, a ,00 + E + Cu-l(s)} (21) 

!/=l/2 

showing the symmetry under the transformation a — > —a. So it is sufficient to consider < a < 1/2. Eq. ( [j"l|) shows 
furthermore that the Casimir energy is independent of the integer part of the flux. 

The analysis of the Casimir energy proceeds using the methods described in detail in |^,|33|,[54|,[l6| and we follow 
their procedure. The starting point is the contour integral representation 

where F(ji,x) describes the eigenvalue equation, 

Ffax) = Jl^ix) - J fa) + ^J^J^x), (23) 

and the integration path T encloses all the zeros of F(fj,,x). 

Shifting the contour to the imaginary axis, one obtains the following representation, valid for 1/2 < SR(s — 1/2) < 1, 

<£*>(•) = ^-i r ( a -V 2 )Bin7r( a -l/2) f°° ^ _ z (s - 1/2) ^ ^ (24) 



V5rT(«) 

where Tint (z I fJ-i u) contains the eigenvalue equation transformed to the imaginary axis, 

T int (z; (x,u) = A log (u~ 2 ^F(n, uj) (25) 

with 

F(/i,«) = I^W + I^M + ^7 M ( U )7 M _!(n). (26) 



Adding and subtracting the uniform asymptotic Debye expansion of F(fJ,,u), see pp| , p^ ] for details, one arrives at 

N 

cm = -d ( !S. w oo + ^nt( S ) + E ^ int) ( s )> ( 27 ) 



z=-l 



where Zi nt (s) is the zeta function with the asymptotic terms subtracted. With t = v j \J v 2 + x 2 , it reads 

Zint(s) = (28) 

E y dx (z 2 - z 2 ) s l^T- mt (z;v + a,x) + T int (z;i>- a,x) - A^(i/,x,t, z)j . (29) 
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On the other hand, the asymptotic terms are, 
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The last equation may be seen as the definition of the coefficients b^j) an d they are easily determined by a simple 
computer program. Most of the b^^ coefficients needed are listed in Appendix A of fusf . 

In order to ensure Z- mi (s) is finite one has to subtract terms up to N = 3 at least. To improve the numerical 
convergence we have chosen N = 5 in the analysis presented later on. 

The asymptotic terms have to be analysed further in order to extract the poles of the theory. They are conveniently 
expressed in terms of the functions |3^Jl^] 

oo r ni —s — b 

f(s;a,b,x) = ^ v a 

V = \jl 

which naturally arise after performing the x- integrals. In detail the full list of useful results is 
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where 

These representations are very suitable for a numerical analysis of the Casimir energy as a function of the mass and 
the flux. 

For small masses, including vanishing mass, as well as for reading off the poles, the expansions in powers of z are 
more suitable. Again, the techniques to obtain these have been described in JhJ . In the present context we find 

.(-t), .R 2 *' 1 ^, ,. n T(n+8-l) 1 /„ , 0n 1 
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n=0 

which makes the pole structure very explicit. 

Before studying this, let us indicate the treatment of the 'extra' contribution in Eq. (^i"|), — Ci/^Iiq,! ( s )- Given no 
/y-sum is involved, it is actually sufficient to just subtract the large argument expansion of Ti nt (z; 1/2 — |a|, x). The 
leading terms in this expansion are 

T int (z; 1/2 - \a\,x) ~ 2 + -LL + + + { X -% (39) 

X x A x 6 
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and a suitable representation to find the analytical continuation in the critical subspace is 
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(40) 



While the first two terms are suitable for a numerical evaluation at s = —1/2, the asymptotic terms are standard 
representations of hypergeometric functions ]3j| ■ 

We are now fully prepared to analyse the pole structure and, if a physically senseful interpretation of finite parts 
is possible, also for a numerical analysis of the Casimir energy. As mentioned, the poles are best read off in the 
representations (|36|) — (j38[) and fl40|). The residues of the single asymptotic terms are 
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the extra contribution adds 



Res I 



c8£,>) 



27Ti? 2 



Summing up, the total residue reads 
Res| s __i Cint(s) = 
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So the pole in the Casimir energy depends on the mass, as to be expected, as well as on the flux via a. This last 
dependence on a did not occur in 2 + 1 dimensions, so that the difference of two situations with different fluxes gave 
a finite answer jl6| . Here this happens only in the case z = so that we should restrict to this case in order to have 
unambiguous answers. In the case z/Oa flux-dependent counterterm FR has to be introduced and finite ambiguities 
remain. The numerical results for the Casimir energy difference between the cases with and without flux for z — are 
presented in Fig. [j] in a black line. Other results, affected by the regularization procedure, are presented as the grey 
lines, corresponding, from top to bottom, to the values of z = 1/64, 1/32, 1/16, 1/8,3/16. For this numerical analysis 
a minimal subtraction was considered. 
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FIG. 1. (Differences of)Interior Casimir Energies - Type I 



Our analysis can be straightforwardly extended to the Type II case. It is not difficult to show that just the 
replacement of |a| by a in the extra term of Eq. (pl|), — C\/2- a { s )% or m tne second order of the asymptotic expansion, 
leads to the wanted results. A similar replacement rule gives the total residue and all the considerations above apply 
directly. The corresponding numerical results are shown in Fig. 
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FIG. 2. (Differences of) Interior Casimir Energies - Type II 



IV. THE CASIMIR ENERGY OF THE EXTERNAL SECTOR 

Let us now consider the external sector in order to see if the pole structure considerably simplifies when combining 
the internal and external spaces. Partly this is to be expected for simple geometric reasons, because the extrinsic 
curvature has the opposite sign. 

The external sector is analysed most effectively using the formulation in terms of a Jost function, see e.g. M. As 
a first step, a large 'external' sphere is introduced as an additional boundary to compactify the space. Boundary 
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conditions are imposed, e.g. MIT ones, and the eigenfunctions are combinations of the ones given in (|4J) and (||) with 
the J^, and replaced by Hu> and . Having the eigenfunctions at hand, it is not difficult to find the Jost 
function for the associated scattering problem. As observed already in various examples |3^ , [l^] , the external sector 
is obtained when I v is replaced by K v . This, essentially, is true also in the present situation apart from the fact that 
the external space does not have a critical subspacc because the singular flux is not sensed directly and no self-adjoint 
extension is involved. The last comment, as it turned out, is just equivalent to the use of a instead of |a| where 
applicable. 

As a result of the above comments, the zeta function of the external part reads 



v=l/2 



where 
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In clear analogy the equations corresponding to ( p7j ) — ( pl| ) can be written down, where one has the relations 



A (ext) _ A (int) 
^-1 — ^-1 ' 



, (ext) 
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and the multipliers are obtained in a similar way as the internal multipliers Also the 'extra' contribution 



"Ci/2-q ^ s dealt with as before, see eq. (|39|) and (f40|), but now we have 



T ext (z; 1/2 - a,x) 



-2-* 
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a z + z 
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The residues in the exterior space are thus obtained from the interior ones by the simple replacement rules already 
mentioned, i.e. z — > — z and \a\ — > a. As a result, the Casimir energy difference between the case with and without 
flux for massless fields is finite and is shown as the black line in Fig. |^. The grey lines show the regularization 
dependent results for, from bottom to top, z = 1/64, 1/32, 1/16, 1/8,3/16. 
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FIG. 3. (Differences of) Exterior Casimir Energies 
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V. THE CASIMIR ENERGY IN THE WHOLE SPACE - 



FINAL COMMENTS 



Adding up the interior and the external sector, we find the pole 
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(50) 



for behaviour I, still showing a flux dependent term. Numerical results for the difference between arbitrary and integer 
flux (only unambiguous for the black z = curve) are shown in Fig. 0. 
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FIG. 4. (Differences of) Complete Casimir Energies - Type I 



On the other hand, 



Res| s ^_i Cint(s) + Res| s ^_i Ccxt(s) = 
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(51) 



gives the residue for behaviour II. We see that the Casimir energy difference between two flux-situations is finite and 
we can fully analyse its dependence on the mass. Numerical results are displayed in Fig. |5|, for the values of z already 
given, from top to bottom. Whereas for non-vanishing mass the field free situation leads to the lowest Casimir energy, 
for non-zero mass a non-trivial global minimum of the Casimir energy emerges. Above a certain critical flux, quantum 
contributions support the stability of the flux, at least for type II. 
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FIG. 5. (Difference of) Casimir Energies - Type II 

In summary, we have studied in detail the behaviour of the Casimir energy of Dirac fields in cylindrically symmetric 
situations. 

In contrast to the 2 + 1-case, when studying the interior case, flux-dependent residues do not vanish when just 
considering the difference between two values of the flux (except for the massless case). In our hope of obtaining 
finite results when considering the whole space, the external space was analysed too. The geometrically grounded 
expectation of a simpler pole structure for the whole space eventually proved to be right, at least for the Type II 
choice of the self-adjoint extension. Completely well-defined results for the massive case were obtained numerically. 

The presence of the mass and the choice of the self-adjoint extension for the radial Hamiltonian at the origin are 
crucial in our results which emphasizes that, for curved boundaries, the contribution of the mass is not exponentially 
small [^3| as it is for parallel plates (37). Concerning the self-adjoint extensions, the remarkable differences in the 
resulting Casimir energies reinforce the well-known fact that the one-parameter family of self-adjoint extensions 
describes nontrivial physics in the core, see also ^2j. Furthermore, a better understanding of the physical meaning 
of the different self-adjoint extensions, by considering how the Casimir energy depends on the parameter of the 
one-parameter family of self-adjoint extensions, should be envisaged in the future. 
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